Grand canonical ensemble approach to electrochemical thermodynamics,kinetics, and model Hamiltonians

The unique feature of electrochemistry is the ability to control reaction thermodynamics and kinetics by the application of electrode potential. Recently, theoretical methods and computational approaches within the grand canonical ensemble (GCE) have enabled to explicitly include and control the electrode potential in first principles calculations. In this review, recent advances and future promises of GCE density functional theory and rate theory are discussed. Particular focus is devoted to considering how the GCE methods either by themselves or combined with model Hamiltonians can be used to address intricate phenomena such as solvent/electrolyte effects and nuclear quantum effects to provide a detailed understanding of electrochemical reactions and interfaces.     

Interplay of Hückel and Möbius Aromaticity in Metal-Metal Quintuple Bonded Complexes of Cr,Mo, and W with Amidinate Ligand: Ab initio DFT and Multireference Analysis**

The aromaticity of metal-metal quintuple bonded complexes of the type M₂L₂ (M=Cr, Mo, and W; L=amidinate) are studied employing gauge including magnetically induced ring current (GIMIC) analysis and electron density of delocalized bonds (EDDB). It is found that the complexes possess two types of aromaticity: i) Hückel aromaticity through delocalization of ligand π electrons with metal-metal δ-bond-forming 6 conjugated electrons (4π and 2δ) ring; ii) Craig-Möbius aromaticity through delocalization of π electrons of both the ligands with metal d-orbitals in Craig type orientation forming 10π electrons ring with a double twist. Extended transition state natural orbital chemical valence (ETS-NOCV) and canonical molecular orbital natural chemical shielding (CMO-NCS) analysis confirm the Craig-Möbius type arrangement of the orbitals. Furthermore, the unprecedented Hückel and Möbius type aromaticity is confirmed from the plot of the current pathways using 3D line integral convolution (3D-LIC) plots. The metal-metal bond order also increases down the group as justified from the complete active space self-consistent field (CASSCF) analysis. Due to an increase in the π and δ electron conjugation, both the Hückel and Möbius aromaticity increase down the group.     

Phosphorene Supported Single-Atom Catalysts for CO Oxidation: A Computational Study

Single-atom catalysts (SACs) have attracted extensive attention owing to their high catalytic activity. The development of efficient SACs is crucial for applications in heterogeneous catalysis. In this article, the geometric configuration, electronic structure, stabilitiy and catalytic performance of phosphorene (Pn) supported single metal atoms (M=Ru, Rh, Pd, Ir, Pt, and Au) have been systematically investigated using density functional theory calculations and ab initio molecular dynamics simulations. The single atoms are found to occupy the hollow site of phosphorene. Among the catalysts studied, Ru-decorated phosphorene is determined to be a potential catalyst by evaluating adsorption energies of gaseous molecules. Various mechanisms including the Eley-Rideal (ER), Langmuir-Hinshelwood (LH) and trimolecular Eley-Rideal (TER) mechanisms are considered to validate the most favourable reaction pathway. Our results reveal that Ru−Pn exhibits outstanding catalytic activity toward CO oxidation reaction via TER mechanism with the corresponding rate-determining energy barrier of 0.44 eV, making it a very promising SAC for CO oxidation under mild conditions. Overall, this work may provide a new avenue for the design and fabrication of two-dimensional materials supported SACs for low-temperature CO oxidation.     

Hyperpolarized 13C Magnetic Resonance Imaging of Fumarate Metabolism by Parahydrogen-induced Polarization: A Proof-of-Concept in vivo Study

Hyperpolarized [1-¹³C]fumarate is a promising magnetic resonance imaging (MRI) biomarker for cellular necrosis, which plays an important role in various disease and cancerous pathological processes. To demonstrate the feasibility of MRI of [1-¹³C]fumarate metabolism using parahydrogen-induced polarization (PHIP), a low-cost alternative to dissolution dynamic nuclear polarization (dDNP), a cost-effective and high-yield synthetic pathway of hydrogenation precursor [1-¹³C]acetylenedicarboxylate (ADC) was developed. The trans-selectivity of the hydrogenation reaction of ADC using a ruthenium-based catalyst was elucidated employing density functional theory (DFT) simulations. A simple PHIP set-up was used to generate hyperpolarized [1-¹³C]fumarate at sufficient ¹³C polarization for ex vivo detection of hyperpolarized ¹³C malate metabolized from fumarate in murine liver tissue homogenates, and in vivo ¹³C MR spectroscopy and imaging in a murine model of acetaminophen-induced hepatitis.     

Dual-Metal Electrolytes for Hybrid-Ion Batteries: Synergism or Antagonism?

The construction of hybrid metal-ion batteries faces a plethora of challenges. A critical one is to unveil the solvation/desolvation processes at the molecular level in electrolytes that ensure efficient transfer of several types of charge carriers. This study reports first results on simulations of mixed-ion electrolytes. All combinations of homo- and hetero-binuclear complexes of Li⁺, Na⁺ and Mg²⁺, solvated with varying number of ethylene carbonate (EC) molecules are modeled in non-polar and polar environment by means of first principles calculations and compared to the mononuclear analogues in terms of stability, spatial organization, charge distribution and solvation/desolvation behavior. The used PF₆⁻ counterion is shown to have minor impact on the geometry of the complexes. The desolvation energy penalty of binuclear complexes can be lowered by the fluoride ions, emerging upon the PF₆⁻ decay. These model investigations could be extended to rationalize the solvation structure and ionic mobility in dual-ion electrolytes.     

Understanding the Chloride Affinity of Barbiturates for Anion Receptor Design

Due to their potential binding sites, barbituric acid (BA) and its derivatives have been used in metal coordination chemistry. Yet their abilities to recognize anions remain unexplored. In this work, we were able to identify four structural features of barbiturates that are responsible for a certain anion affinity. The set of coordination interactions can be finely tuned with covalent decorations at the methylene group. DFT-D computations at the BLYP-D3(BJ)/aug-cc-pVDZ level of theory show that the C−H bond is as effective as the N−H bond to coordinate chloride. An analysis of the electron charge density at the C−H⋅⋅⋅Cl⁻ and N−H⋅⋅⋅Cl⁻ bond critical points elucidates their similarities in covalent character. Our results reveal that the special acidity of the C−H bond shows up when the methylene group moves out of the ring plane and it is mainly governed by the orbital interaction energy. The amide and carboxyl groups are the best choices to coordinate the ion when they act together with the C−H bond. We finally show how can we use this information to rationally improve the recognition capability of a small cage-like complex that is able to coordinate NaCl.     

Rational Development of a Metal-Free Bifunctional System for the C−H Activation of Methane: A Density Functional Theory Investigation

The activation or heterolytic splitting of methane, a challenging substrate usually restricted to transition metals, has so far proven elusive in experimental frustrated Lewis pair (FLP) chemistry. In this article, we demonstrate, using density functional theory (DFT), that 1-aza-9-boratriptycene is a conceptually simple intramolecular FLP for the activation of methane. Systematic comparison with other FLP systems allows to gain insight into their reactivity with methane. The thermodynamics and kinetics of methane activation are interpreted by referring to the analysis of the natural charges and by employing the distortion-interaction/activation strain (DIAS) model. These showed that the nature of the Lewis base influences the selectivity over the reaction pathway, with N Lewis bases favoring the deprotonation mechanism and P bases the hydride abstraction one. The lower barrier of activation for 1-aza-9-boratriptycene and the higher products stability are due to a better interaction energy than its counterparts, itself due to electrostatic interactions with the methane moiety, favorable orbital overlaps allowed by the side-attack, and space proximity between the B and N atoms.     

Functional Group Modification and Bonding Characteristics of Ti3C2 MXene-Organic Composites from First-Principles Calculations

The unique physical structure and abundant surface functional groups of MXene make the grafted organic molecules exhibit specific electrical and optical properties. This work reports the results of first-principles calculations to investigate the composite systems formed by different organic molecular monomers, namely acrylic acid (AA), acrylamide (AM), 1-aziridineethanol (1-AD) and glucose, and Ti₃C₂ MXene saturated with different functional groups, namely −OH, −O and −F. The results show that the interaction between organic molecules and the MXene surface depends on the type of functional groups of the organic molecules, while the strength of the interaction is determined by the type of surface functional groups and the number of hydrogen bonds. The bare Ti₃C₂ and Ti₃C₂(OH)₂ can readily form strong chemical and hydrogen bonds with AA and AM molecules, leading to strong adsorption energy and a large amount of charge transfer, while the interaction between organic molecules and MXene saturated by −F or −O groups mainly exhibits physical interactions, accompanied by low adsorption energy and a small amount of charge transfer. This research provides theoretical guidance for the synthesis of high-performance MXene organic composite systems.     

Indirectly Detected DNP-Enhanced 17O NMR Spectroscopy: Observation of Non-Protonated Near-Surface Oxygen at Naturally Abundant Silica and Silica-Alumina

Recent studies have shown that dynamic nuclear polarization (DNP) can be used to detect ¹⁷O solid-state NMR spectra of naturally abundant samples within a reasonable experimental time. Observations using indirect DNP, which relies on ¹H mediation in transferring electron hyperpolarization to ¹⁷O, are currently limited mostly to hydroxyls. Direct DNP schemes can hyperpolarize non-protonated oxygen near the radicals; however, they generally offer much lower signal enhancements. In this study, we demonstrate the detection of signals from non-protonated ¹⁷O in materials containing silicon. The sensitivity boost that made the experiment possible originates from three sources: indirect DNP excitation of ²⁹Si via protons, indirect detection of ¹⁷O through ²⁹Si nuclei using two-dimensional ²⁹Si{¹⁷O} D-HMQC, and Carr-Purcell-Meiboom-Gill refocusing of ²⁹Si magnetization during acquisition. This ²⁹Si-detected scheme enabled, for the first time, 2D ¹⁷O−²⁹Si heteronuclear correlation spectroscopy in mesoporous silica and silica-alumina surfaces at natural abundance. In contrast to the silanols showing motion-averaged ¹⁷O signals, the framework oxygens exhibit unperturbed powder patterns as unambiguous fingerprints of surface sites. Along with hydroxyl oxygens, detection of these moieties will help in gaining more atomistic-scale insights into surface chemistry.     

Computational Screening of 3 d Transition Metal Atoms Anchored on Defective Graphene for Efficient Electrocatalytic N2 Fixation

The synthesis of ammonia (NH₃) through the electrochemical reduction of molecular nitrogen (N₂) is a promising strategy for significantly reducing energy consumption compared to traditional industrial processes. Herein, we report the design of a series of monovacancy and divacancy defective graphenes decorated with single 3d transition metal atoms (TM@MVG and TM@DVG; TM=Sc−Zn) as electrocatalysts for the nitrogen-reduction reaction (NRR) aided by density functional theory (DFT) calculations. By comparing energies for N₂ adsorption as well as the free energies associated with *N₂ activation and *N₂H formation, we successfully identified V@MVG, with the lowest potential of −0.63 V, to be an effective catalytic substrate for the NRR in an enzymatic mechanism. Electronic properties, including Bader charges, charge density differences, partial densities of states, and crystal orbital Hamilton populations, are further analyzed in detail. We believe that these results help to explain recent observations in this field and provide guidance for the exploration of efficient electrocatalysts for the NRR.